Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression.

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

An immunogenetic basis for lung cancer risk.

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants.

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Validation of a Machine Learning Model to Predict Immunotherapy Response in Head and Neck Squamous Cell Carcinoma.

Head and neck squamous-cell carcinoma (HNSCC) is a disease with a generally poor prognosis; half of treated patients eventually develop recurrent and/or metastatic (R/M) disease. Patients with R/M HNSCC generally have incurable disease with a median survival of 10 to 15 months. Although immune-checkpoint blockade (ICB) has improved outcomes in patients with R/M HNSCC, identifying patients who are likely to benefit from ICB remains a challenge. Biomarkers in current clinical use include tumor mutational burden and immunohistochemistry for programmed death-ligand 1, both of which have only modest predictive power. Machine learning (ML) has the potential to aid in clinical decision-making as an approach to estimate a tumor's likelihood of response or a patient's likelihood of experiencing clinical benefit from therapies such as ICB. Previously, we described a random forest ML model that had value in predicting ICB response using 11 or 16 clinical, laboratory, and genomic features in a pan-cancer development cohort. However, its applicability to certain cancer types, such as HNSCC, has been unknown, due to a lack of cancer-type-specific validation. Here, we present the first validation of a random forest ML tool to predict the likelihood of ICB response in patients with R/M HNSCC. The tool had adequate predictive power for tumor response (area under the receiver operating characteristic curve = 0.65) and was able to stratify patients by overall (HR = 0.53 [95% CI 0.29-0.99], p = 0.045) and progression-free (HR = 0.49 [95% CI 0.27-0.87], p = 0.016) survival. The overall accuracy was 0.72. Our study validates an ML predictor in HNSCC, demonstrating promising performance in a novel cohort of patients. Further studies are needed to validate the generalizability of this algorithm in larger patient samples from additional multi-institutional contexts.

Pre-treatment serum albumin and mutational burden as biomarkers of response to immune checkpoint blockade.

The effects of cytokine and protein stabilizing carriers, such as serum albumin, on tumor response to immune checkpoint blockade (ICB) is not well understood. By examining 1714 patients across 16 cancer types, we found that high pretreatment serum albumin level predicts favorable tumor radiographic response following ICB treatment in a dose-dependent fashion. Serum albumin is a candidate biomarker that can be combined with tumor mutational burden (TMB) for additional predictive capacity, and the tumor response rate to ICB was ~49% in the albumin-high/TMB-high group.

Analysis of surface characteristics of (Y, Nb)-TZP after finishing and polishing.

PURPOSE: This in vitro study aimed to evaluate the surface characteristics of a full veneer crown fabricated chairside (CS) from a (Y, Nb)-TZP zirconia block in response to conventional zirconia grinding and polishing. MATERIALS AND METHODS: Zirconia crowns (n = 40) were first prepared and divided into two groups of materials: Labside (LS) and CS, after which each specimen went through a five-step grinding and polishing procedure. Following each surface treatment, surface characteristics were analyzed using confocal laser microscopy (CLSM), average surface roughness (Ra) values were processed from the profile data through Gaussian filtering, and X-ray diffraction pattern analysis was performed to evaluate the monoclinic (M) phase content. Then, a representative specimen was selected for field-emission scanning electron microscopy (FE-SEM), followed by a final analysis of the roughness and X-ray diffraction of the specimens using the independent t-test and repeated measures analysis of variance (RM-ANOVA). RESULTS: In every group, polishing significantly reduced the Ra values (P < .001). There was no significant difference in Ra between the polished state CS and LS. Furthermore, CLSM and FE-SEM investigations revealed that even though grain exposure was visible in CS specimens throughout the as-delivered and ground states, the exposure was reduced after polishing. Moreover, while no phase transformation was visible in the LS, phase transformation was visible in CS after every surface treatment, with the M phase content of the CS group showing a significant reduction after polishing (P < .001). CONCLUSION: Within the limits of this study, clinically acceptable level of surface finishing of (Y, Nb)-TZP can be achieved after conventional zirconia polishing sequence.

Improved prediction of immune checkpoint blockade efficacy across multiple cancer types.

Only a fraction of patients with cancer respond to immune checkpoint blockade (ICB) treatment, but current decision-making procedures have limited accuracy. In this study, we developed a machine learning model to predict ICB response by integrating genomic, molecular, demographic and clinical data from a comprehensively curated cohort (MSK-IMPACT) with 1,479 patients treated with ICB across 16 different cancer types. In a retrospective analysis, the model achieved high sensitivity and specificity in predicting clinical response to immunotherapy and predicted both overall survival and progression-free survival in the test data across different cancer types. Our model significantly outperformed predictions based on tumor mutational burden, which was recently approved by the U.S. Food and Drug Administration for this purpose1. Additionally, the model provides quantitative assessments of the model features that are most salient for the predictions. We anticipate that this approach will substantially improve clinical decision-making in immunotherapy and inform future interventions.

Indigenous Ancestry and Admixture in the Uruguayan Population.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in "Charrúas" and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

Genomic and Transcriptomic Characteristics According to Size of Papillary Thyroid Microcarcinoma.

It is controversial as to whether papillary thyroid microcarcinoma (PTMC) has some genomic and transcriptomic characteristics that differentiate between an early-stage lesion that would eventually evolve into the larger papillary thyroid cancer (PTC), and an occult indolent cancer in itself. To investigate this, we comprehensively elucidated the genomic and transcriptomic landscapes of PTMCs of different sizes, using a large-scaled database. This study included 3435 PTCs, 1985 of which were PTMCs. We performed targeted next-generation sequencing for 221 PTCs and integrated these data with the data including The Cancer Genome Atlas (TCGA) project. The frequency of v-raf murine sarcoma viral oncogene homolog B (BRAF)V600E mutation was higher in PTMCs >0.5 cm than that in very small PTMCs (≤0.5 cm) and decreased again in PTCs >2 cm. Among PTMCs, the prevalence of mutations in rat sarcoma (RAS) and telomerase reverse transcriptase (TERT) promoter was not significantly different according to their size, but lower than in large PTCs. There was no change in the tumor mutational burden, the number of driver mutations, and transcriptomic profiles with tumor size, among PTMCs and all PTCs. Although a few genes with differential expression and TERT promoter mutations were found in a few PTMCs, our findings showed that there were no useful genomic or transcriptomic characteristics for the prediction of the future progression of PTMC.

Recent Improvements in Genomic and Transcriptomic Understanding of Anaplastic and Poorly Differentiated Thyroid Cancers.

Anaplastic thyroid cancer (ATC) is a lethal human cancer with a 5-year survival rate of less than 10%. Recently, its genomic and transcriptomic characteristics have been extensively elucidated over 5 years owing to advance in high throughput sequencing. These efforts have extended molecular understandings into the progression mechanisms and therapeutic vulnerabilities of aggressive thyroid cancers. In this review, we provide an overview of genomic and transcriptomic alterations in ATC and poorly-differentiated thyroid cancer, which are distinguished from differentiated thyroid cancers. Clinically relevant genomic alterations and deregulated signaling pathways will be able to shed light on more effective prevention and stratified therapeutic interventions for affected patients.

BRAFV600E Transduction of an SV40-Immortalized Normal Human Thyroid Cell Line Induces Dedifferentiated Thyroid Carcinogenesis in a Mouse Xenograft Model.

Background: Despite active studies of the clinical importance of BRAFV600E, suitable research models to investigate the role of this mutation in the etiopathogenesis of human thyroid cancers are limited. Thus, we generated cell lines by transducing the simian virus (SV)-40 immortalized human thyroid cell line Nthy-ori 3-1 (Nthy) with lentiviral vectors expressing either BRAFWT (Nthy/WT) or BRAFV600E. Nthy/WT and Nthy/V600E cells were then xenografted into mice to evaluate the carcinogenic role of BRAFV600E. Methods: Each cell line was subcutaneously injected into NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ mice, and a pathological analysis was performed. The effects of the mutation were further verified by using a BRAFV600E-selective inhibitor (PLX-4032, vemurafenib). The transcriptome was analyzed by RNA sequencing and compared with data from The Cancer Cell Line Encyclopedia and Gene Expression Omnibus. Results: While Nthy/WT was not tumorigenic in vivo, Nthy/V600E formed tumors reaching 2784.343 mm3 in 4 weeks, on average. A pathological analysis indicated that Nthy/V600E tumors were dedifferentiated thyroid cancer. We found metastases in the lung, liver, and relevant lymph nodes. A transcriptomic analysis revealed 5512 differentially expressed genes (DEGs) between the mutant and wild-type cell lines, and more DEGs were shared with anaplastic thyroid cancer than with papillary thyroid cancer. BRAFV600E activated the cell cycle mainly by regulating G1/S phases. PLX-4032 treatment significantly inhibited tumor growth and metastasis. Conclusions: Our data show that BRAFV600E plays a pivotal role in the carcinogenic transformation of an SV40-transfected immortalized normal human thyroid cell line. This xenograft model is expected to contribute to studies of the etiopathogenesis and treatment of highly malignant thyroid cancers.

NARD: whole-genome reference panel of 1779 Northeast Asians improves imputation accuracy of rare and low-frequency variants.

Here, we present the Northeast Asian Reference Database (NARD), including whole-genome sequencing data of 1779 individuals from Korea, Mongolia, Japan, China, and Hong Kong. NARD provides the genetic diversity of Korean (n = 850) and Mongolian (n = 384) ancestries that were not present in the 1000 Genomes Project Phase 3 (1KGP3). We combined and re-phased the genotypes from NARD and 1KGP3 to construct a union set of haplotypes. This approach established a robust imputation reference panel for Northeast Asians, which yields the greatest imputation accuracy of rare and low-frequency variants compared with the existing panels. NARD imputation panel is available at https://nard.macrogen.com/ .

CXCL16 positively correlated with M2-macrophage infiltration, enhanced angiogenesis, and poor prognosis in thyroid cancer.

Although various chemokines have pro-tumorigenic actions in cancers, the effects of CXCL16 remain controversial. The aim of this study was to investigate the molecular characteristics of CXCL16-expressing papillary thyroid cancers (PTCs). CXCL16 expressions were significantly higher in PTCs than benign or normal thyroid tissues. In the TCGA dataset for PTCs, a higher CXCL16 expression was associated with M2 macrophage- and angiogenesis-related genes and poor prognostic factors including a higher TNM staging and the BRAFV600E mutation. PTCs with a higher expression of 3-gene panel including CXCL16, AHNAK2, and THBS2 showed poor recurrence-free survivals than that of the lower expression group. Next, shCXCL16 was introduced into BHP10-3SCp cells to deplete the endogenous CXCL16, and then, the cells were subcutaneously injected to athymic mice. Tumors from the BHP10-3SCpshCXCL16 exhibited a delayed tumor growth with decreased numbers of ERG+ endothelial cells and F4/80+ macrophages than those from the BHP10-3SCpcontrol. CXCL16-related genes including AHNAK2 and THBS2 were downregulated in the tumors from the BHP10-3SCpshCXCL16 compared with that from the BHP10-3SCpcontrol. In conclusion, a higher CXCL16 expression was associated with macrophage- and angiogenesis-related genes and aggressive phenotypes in PTC. Targeting CXCL16 may be a good therapeutic strategy for advanced thyroid cancer.

Integrative analysis of genomic and transcriptomic characteristics associated with progression of aggressive thyroid cancer.

Anaplastic thyroid cancer (ATC) and advanced differentiated thyroid cancers (DTCs) show fatal outcomes, unlike DTCs. Here, we demonstrate mutational landscape of 27 ATCs and 86 advanced DTCs by massively-parallel DNA sequencing, and transcriptome of 13 ATCs and 12 advanced DTCs were profiled by RNA sequencing. TERT, AKT1, PIK3CA, and EIF1AX were frequently co-mutated with driver genes (BRAFV600E and RAS) in advanced DTCs as well as ATC, but tumor suppressors (e.g., TP53 and CDKN2A) were predominantly altered in ATC. CDKN2A loss was significantly associated with poor disease-specific survival in patients with ATC or advanced DTCs, and up-regulation of CD274 (PD-L1) and PDCD1LG2 (PD-L2). Transcriptome analysis revealed a fourth molecular subtype of thyroid cancer (TC), ATC-like, which hardly reflects the molecular signatures in DTC. Furthermore, the activation of JAK-STAT signaling pathway could be a potential druggable target in RAS-positive ATC. Our findings provide insights for precision medicine in patients with advanced TCs.

Interaction of BRAF-induced ETS factors with mutant TERT promoter in papillary thyroid cancer.

Synergistic effects of BRAFV600E and TERT promoter mutations on the poor clinical outcomes in papillary thyroid cancer (PTC) have been demonstrated. The potential mechanism of this phenomenon has been proposed: MAPK pathway activation by the BRAFV600E mutation may upregulate E-twenty six (ETS) transcription factors, increasing TERT expression by binding to the ETS-binding site generated by the TERT promoter mutation; however, it has not yet been fully proven. This article provides transcriptomic insights into the interaction between BRAFV600E and TERT promoter mutations mediated by ETS factors in PTC. RNA sequencing data on 266 PTCs from The Cancer Genome Atlas and 65 PTCs from our institute were analyzed for gene expression changes and related molecular pathways, and the results of transcriptomic analyses were validated by in vitro experiments. TERT mRNA expression was increased by the coexistence of BRAFV600E and TERT promoter mutations (fold change, 16.17; q-value = 7.35 × 10-12 vs no mutation). In the ETS family of transcription factors, ETV1, ETV4 and ETV5 were upregulated by the BRAFV600E/MAPK pathway activation. These BRAFV600E-induced ETS factors selectively bound to the mutant TERT promoter. The molecular pathways activated by BRAFV600E were further augmented by adding the TERT promoter mutation, and the pathways related to immune responses or adhesion molecules were upregulated by TERT expression. The mechanism of the synergistic effect between BRAFV600E and TERT promoter mutations on cancer invasiveness and progression in PTC may be explained by increased TERT expression, which may result from the BRAF-induced upregulation of several ETS transcription factors.

Comprehensive Transcriptomic and Genomic Profiling of Subtypes of Follicular Variant of Papillary Thyroid Carcinoma.

BACKGROUND: Among subtypes of follicular variant of papillary thyroid carcinoma (FVPTC), encapsulated FVPTC (EFVPTC) shows more indolent behavior than infiltrative FVPTC (IFVPTC). In particular, noninvasive EFVPTC, now designated as noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP), tends to have an excellent prognosis. However, it remains unclear whether the molecular pathogenesis or signature of the various forms of FVPTC is different. By massively parallel sequencing analysis, this study comprehensively characterized the transcriptional and mutational landscape of FVPTC and established correlations with phenotypic subtypes. METHODS: This study included 48 FVPTCs: 17 NIFTPs, 13 invasive EFVPTCs (I-EFVPTCs), and 18 IFVPTCs. For comparison, 55 classical papillary thyroid carcinomas (cPTCs) harboring a BRAFV600E mutation, six follicular adenomas (FAs), and 15 minimally invasive follicular thyroid carcinomas (miFTCs) with RAS mutations were also included. RESULTS: In NIFTP, the BRAFV600E mutation was not found, but RAS and other alterations were present in 64.7% and 17.6% of cases, respectively. However, in I-EFVPTC and IFVPTC, the proportions of BRAFV600E mutation (38.5% and 38.9%, respectively) and of RAS mutations (38.5% and 38.9%, respectively) or other alterations (15.4% and 16.7%, respectively) were similar. On a molecular level, RAS-mutated FVPTCs were all RAS-like except for one IFVPTC case. Transcriptomic profiles of NIFTP, I-EFVPTC, and FA/miFTC were comparable, although the profile of RAS-mutated IFVPTC was altered to activate molecular pathways involved in cell adhesion and invasion. Interestingly, 80% of BRAFV600E-mutated I-EFVPTCs were also classified as RAS-like, whereas all BRAFV600E-mutated IFVPTCs were BRAF-like and indistinguishable from cPTC. Molecular pathways associated with cell adhesion and invasion were also differentially activated in BRAFV600E-mutated IFVPTC. CONCLUSIONS: Molecular profiles of NIFTP and I-EFVPTC may be shared with FA/miFTC, while IFVPTC seems to be associated with a similar profile as cPTC. Activation of cell adhesion and invasion pathways may play a key role in the development of invasive phenotypes of FVPTC.

Signatures of photo-aging and intrinsic aging in skin were revealed by transcriptome network analysis.

There are various factors that alter physiological characteristics in skin. Elucidating the underlying mechanism of transcriptional alterations by intrinsic and extrinsic factors may lead us to understand the aging process of skin. To identify the transcriptomic changes of the aging skin, we analyzed publicly available RNA sequencing data from Genotype-Tissue Expression (GTEx) project. GTEx provided RNA sequencing data of suprapubic (n=228) and lower leg (n=349) skins, which are photo-protected and photo-damaged. Using differentially expressed gene analysis and weighted gene co-expression network analysis, we characterized transcriptomic changes due to UV exposure and aging. Genes involved in skin development such as epidermal differentiation complex component (SPRR and LCE families), vasculature development (TGFBR1, TGFBR2, TGFBR3, KDR, FGF2, and VEGFC), and matrix metalloproteinase (MMP2, MMP3, MMP8, MMP10, and MMP13) were up-regulated by UV exposure. Also, down-regulated lipid metabolism and mitochondrial biogenesis were observed in photo-damaged skin. Moreover, wound healing process was universally down-regulated in suprapubic and lower leg with aging and further down-regulation of lipid metabolism and up-regulation of vasculature development were found as photo-aging signatures. In this study, dynamic transcriptomic alterations were observed in aged skin. Hence, our findings may help to discover a potential therapeutic target for skin rejuvenation.

Development of a common platform for the noninvasive prenatal diagnosis of X-linked diseases.

OBJECTIVE: The aim of this study was to develop a common targeted massively parallel sequencing platform for the noninvasive prenatal diagnosis (NIPD) of multiple X-linked diseases. METHOD: The custom capture probe was designed to target 33 genes and recombination hotspots. We tested the carrier mother and male proband pair of 6 families. Plasma DNA of the pregnant carrier mother was collected at different gestational weeks and sequenced. The fetal genotype of each family was determined by estimating the imbalance between the 2 maternal haplotypes constructed using a common custom-designed platform. RESULTS: The targeted sequencing of the maternal, proband, and fetal genomic DNAs and maternal plasma DNAs resulted in uniform coverage across the target region. Three to 5 recombination points were observed in each sample. However, these recombination points did not affect the haplotype dosage analysis for fetal genotype prediction. Consequently, all fetal genotypes in the 6 families obtained from haplotype dosage analysis of maternal plasma sequencing data were predicted correctly. CONCLUSIONS: Since a single platform that covers multiple diseases may prevent the need for disease-specific probes for the NIPD of individual disorders, this approach may provide a practical advantage for clinically implementing the NIPD of X-linked diseases.

Targeted linked-read sequencing for direct haplotype phasing of maternal DMD alleles: a practical and reliable method for noninvasive prenatal diagnosis.

For the noninvasive prenatal diagnosis (NIPD) of X-linked recessive diseases such as Duchenne muscular dystrophy (DMD), maternal haplotype phasing is a critical step for dosage analysis of the inherited allele. Until recently, the proband-based indirect haplotyping method has been preferred despite its limitations for use in clinical practice. Here, we describe a method for directly determining the maternal haplotype without requiring the proband's DNA in DMD families. We used targeted linked-read deep sequencing (mean coverage of 692×) of gDNA from 5 mothers to resolve their haplotypes and predict the mutation status of the fetus. The haplotype of DMD alleles in the carrier mother was successfully phased through a targeted linked-read sequencing platform. Compared with the proband-based phasing method, linked-read sequencing was more accurate in differentiating whether the recombination events occurred in the proband or in the fetus. The predicted inheritance of the DMD mutation was diagnosed correctly in all 5 families in which the mutation had been confirmed using amniocentesis or chorionic villus sampling. Direct haplotyping by this targeted linked-read sequencing method could be used as a phasing method for the NIPD of DMD, especially when the genomic DNA of the proband is unavailable.